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In this contribution we demonstrate structural control over a transport resonance in HS(CH 2 )n [1,4 
- C 6 H 4 ](CH 2 )nSH (n = 1, 3, 4, 6) metal |molecule |metal junctions, fabricated and tested using the 
scanning tunnelling microscopy-based I[z) method. The Breit-Wigner resonance originates from 
one of the arene 7r-bonding orbitals, which sharpens and moves closer to the contact Fermi energy 
as n increases. Varying the number of methylene groups thus leads to a very shallow decay of the 
conductance with the length of the molecule. We demonstrate that the electrical behaviour observed 
here can be straightforwardly rationalized by analyzing the effects caused by the electrostatic bal¬ 
ance created at the metal-molecule interface. Such resonances offer future prospects in molecular 
electronics in terms of controlling charge transport over longer distances, and also in single molecule 
conductance switching if the resonances can be externally gated. 

PACS numbers: 73.63.Rt, 73.40.Sx, 85.65.+h 

I. INTRODUCTION 


The development of reliable techniques for the fab¬ 
rication and electrical testing of metal | single molecule 
I metal junctions, supported by rapid developments in 
theoretical treatment of such junctions, has led to a 
significant re-iuA^oration of the wider field of molecu¬ 
lar electronics.^^^ The techniques used can be divided 
into two categories. One category uses solid-state elec¬ 
tronics technology to fabricate so-called break junctions, 
in which a very thin metallic wire is deliberately bro¬ 
ken (either by electr emigration,!^ or by the application 
of mechanical stres^ in such a way that a molecule 
can bridge the resulting gap. The other category in¬ 
volves use of scanning probe microscopy to fabricate and 
characterise metal |molecule |metal junctions. For exam¬ 
ple, in the scanning tunnelling microscopy break junc¬ 
tion (STM-BJ) method, a gold STM tip is driven into 
a gold substrate, and then retracted while the current 
is monitored.!^ As the resulting chains of gold atoms are 
broken, step decreases in the current approximately cor¬ 
responding to the quantum of conductance, Go, are seen 
until the final chain of gold atoms is broken. In the pres¬ 
ence of molecules that contain terminal po tentia l binding 
groups (e.g. thio ls,!^ pyridines^^ aminescarboxy- 

late anions, phosphine^^^^ , thioetherJ^, addi¬ 
tional much smaller current step decreases are seen after 
the final gold atomic chain breaks, and these correspond 
to the subsequent breaking of metal | molecule | metal 
junctions. Alternatively, in the so-called I{z) (/, current; 
2 :, vertical height) method for measuring single molecule 
electricla properties, a gold STM tip is brought into close 
proximity of a surface coated with suitable molecules 
without allowing it to come into contact, and is then 
retracted while the tunnelling current is monitored.!^ 


(a) 



FIG. 1. Structural comparison of (a) 1,12-dodecanedithiol 
with (b) 6Ph6 and its derivatives (X=2,3,5,6-tetrafluoro-,2,5- 
dimethyl- 2,5-dimet hoxy-) 


When a molecule bridges the gap prior to tip retraction, 
a characteristic plateau is seen in the current-distance 
curve, prior to junction breakdown whereupon the cur¬ 
rent drops sharply. 

For the great majority of (relatively short) molecules 
that have been studied in metal | single molecule | metal 
junctions, the mechanism of conductance has been found 
to be coherent tunnelling.!^! The conductance, Gm, then 
varies as a function of the length of the bridging molecule 
according to Eq. 0 : 

Gm = ( 1 ) 

where Gcon =GoTLTi^ (Go is the quantum unit of 
conductance, 77.4 /iS, and T^ and Tr are the trans- 
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mittances of the left hand and right hand connections 
of the molecule to the metal contacts; Gcon is called the 
contact conductance and represents the conductance of 
the system in the absence of the molecular backbone), R 
is the molecular backbone length and /3 is the attenua¬ 
tion factor. A plot of IuGm vs. R is used to determine 
[ 3 . Often, such molecules can be considered as “passive” 
molecular wires, in which the conductance decreases with 
bridge length in the simple manner predicted by phase 
coherent tunnelling. In many cases this length decay 
can be straightforwardly linked to the degree of conju¬ 
gation, in turn related to the proximity of the relevant 
HOMO or LUMO (“frontier”) orbitals to the metal con¬ 
tact Fermi energies However, molecular electronics 
is now aiming to achieve more active control over the 
electrical properties of molecular bridges. Approaches to 
this include electrochemical switching or conformational 
control of bridge structure and electronic states and con¬ 
trol of quantum interference in molecular circuits. 

In earlier work, we showed that the conductance of 
junctions involving 6Ph6 and its derivatives (Fig. ) 
was much higher ( 0.7 nS) than that measured for 1,12- 
dodecanedithiol (i.e. the two alkyl links for 6Ph6 con¬ 
nected back-to-back; 0.028 nS^, even though 6Ph6 
(2.10 nm S...S distance) is considerably longer than 1,12- 
dodecanedithiol ( 1.72 nm S...S). We suggested that the 
reason for the unusually high conductance of 6Ph6 was 
that it behaves as a ‘molecular double tunnelling barrier’ 
in which a ‘well’ with lower HOMO-LUMO (tt-tt*) sepa¬ 
ration is sandwiched between two ‘barriers’ with higher 
HOMO-LUMO (cr-cr*) separation, by analogy with inor¬ 
ganic double tunnelling barrier structures involving HI- 
V semiconductors . 1 ^ Clearly, the benzene moiety is ca¬ 
pable of acting as an indentation in the tunnelling bar¬ 
rier. Interestingly, control over the conductance of the 
6Ph6 family could be exerted chemically by using sub¬ 
stituents on the aryl ring (Fig. [^; electron-donating sub¬ 
stituents (Me—, MeO—) resulted in significantly higher 
junction conductances and the electron-deficient C6F4 
group resulted in a lower conductance.^ Many studies of 
metal |molecule |metal junctions have probed the molec¬ 
ular structural requirements for obtaining high junction 
conduct ance.fl^M^ 

The work described in this contribution is aimed at 
a deeper understanding of the evolution of the orbital 
alignment in single-molecule junctions. To this end, we 
have investigated the effect of varying the width of the 
‘barriers’ (i.e. the number of methylene groups) while 
keeping the ‘well’ (the 1 , 4 —C6H4 moiety) constant. Ac¬ 
cordingly, we have now synthesised new molecules and 
have determined the single molecule conductances of 
1,4-HS(CH2)^C6H4(CH2)^SH (xPhx; x = 1 , 3 , 4 , 6) us¬ 
ing the I{z) technique, and we have carried out transport 
calculations on the molecules. We find that the variation 
of molecular conductance with n in this family is sur¬ 
prisingly small, and that this is the consequence of the 
presence of a Breit-Wigner resonance in the transmis¬ 
sion curve (caused by the higher-energy of the two tt - 


bonding orbitals of the benzene ring). This resonance 
sharpens and moves closer to the contact Fermi energy 
as n increases, offsetting to some degree the rapid expo¬ 
nential decrease in conductance with length that might 
be expected for increasing n. Thus, control over the 
conductance of these molecules can be exerted not only 
by substituents on the phenyl ring, but also by vary¬ 
ing the alkyl chain length. More importantly, we show 
that this system can be understood in terms of a sim¬ 
ple pair of back-to-back metal-semiconductor junctions 
where the alignment of the states in the internal part of 
the semiconductor primarily depends on how the electro¬ 
static balance is established at the contact region. Such a 
balance is in turn mainly determined between the molec¬ 
ular mid-gap states and the metal Fermi level via charge 
transfer, which finally determines the position of the bulk 
resonance with respect to the Fermi level. We also show 
that investigating this system within this perspective is 
essential to understand the level alignment in the junc¬ 
tion, which otherwise could not be fully explained by the 
initial trend in the gas phase. 

II. EXPERIMENTAL RESULTS 

A. Syntheses, sample preparation and monolayer 
formation 

IPhl (as the free di-thiol) is commercially-available 
and was used as received. 6Ph6 (dithioacetate; Fig. 
was synthesised as previously describee?^ and the same 
general approach was used to synthesise 3 Ph 3 and 4 Ph 4 
; full details are provided in the Supplemental Mate¬ 
rial (SM).I^ Gold-coated glass substrates (Arrandee) of 
ca. 1x1 cm size were flame-annealed using a butane- 
air torch to red heat, to generate Au(lll) terraces,!^ 
as confirmed by STM imaging (see SM, 1 ^ Fig.l). The 
slides were allowed to cool, and were then immersed in a 
10 “^ M solution of the appropriate di-thioacetate ( 3 Ph 3 , 
4 Ph 4 , 6Ph6) or dithiol (IPhl) in GH2GI2 for 1-2 min¬ 
utes (IPhl; 3 minutes) to afford a partial monolayer cov¬ 
erage of adsorbed dithiol molecules. Polarisation mod¬ 
ulated infrared reflection-absorption spectroscopy (PM- 
IRRAS) was used to confirm that the adsorption of the 
molecules took place (for an example, see Fig .2 in the 
SM), and STM imaging of these submonolayers showed 
bright spots indicative of the presence of molecule(s), 
particularly at step edges (see SM,I^ Fig.lc); molecular 
resolution could not be obtained, as is typical for highly- 
mobile adsorbed thiols at submonolayer coverages. 


B. Single molecule conductance determination 

The scanning tunnelling microscopy-based I{z) tech¬ 
nique was used to determine single molecule conductance 
values for IPhl, 3 Ph 3 , 4 Ph 4 and 6Ph6. In this tech¬ 
nique, tip-substrate contact is avoided. A freshly-cut An 
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Conductance (nS) Retraction-► 


FIG. 2. (a) Histogram of 500 I(s) scans obtained for 4Ph4 
using a set point current of 10 nA, and (b) a 2-dimensional 
histogram plot, colour-coded to show frequency of occurrence 
of conductance as a function of distance in these I{s) traces. 



FIG. 3. (a) Histogram of 500 I{s) scans obtained for IPhl 
using a set point current of 30 nA, and (b) a 2-dimensional 
histogram plot, colour-coded to show frequency of occurrence 
of conductance as a function of distance in these I(s) traces. 


STM tip was brought into close proximity to the modi¬ 
fied substrate surface using the appropriate set point cur¬ 
rent, and the tip was then retracted while the tunnelling 
current was measured. When a molecule (or molecules) 
bridge between tip and substrate, characteristic current 
plateaux in the current-vertical distance {I — z) traces re¬ 
sult, and the molecular conductance is determined from a 
histogram of many such traces. In this work, at least 500 
I — z traces showing current plateaux (defined as hav¬ 
ing length > 0.5 nm) were combined into conductance 
histograms for any given experiment; typically, approxi¬ 
mately 15 % of / — 2 : traces showed a plateau. Of the re¬ 
maining traces, approximately 60% were straightforward 
exponential decays with no evidence of molecular bridge 
formation and the rest were noisy and showed no obvi¬ 
ous plateau; I — z traces in these categories were rejected. 
Two-dimensional histograms of conductance-vertical dis¬ 
tance correlation were also constructed. The tip bias in 
all experiments was +0.6 V. 

Initially, we began by repeating the experiment of 
Leary et al.^ using 6Ph6, with a set point current of 
6 nA. This produced a relatively sharp histogram peak 
at a conductance of 0.82 ± 0.18 nS, in good agreement 
with the results reported by Leary et al (see SM+^Fig.4). 
A second clear peak can also be seen at approximately 
twice this value, corresponding to junctions formed with 


Molecule 

Gonductance 

(nS) 

95th percentile 
break-off 
distance (nm) 

Galcd. An...An 
distance for 
An xPhx An 
junction (nm)“ 

IPhl 

4.14 ± 0.58 

1.4 

1.1 

3Ph3 

2.01 ± 0.32 

1.8 

1.7 

4Ph4 

1.21 ± 0.19 

1.9 

1.9 

6 Ph6 

0.82 ± 0.18 

2.5 

2.3 


TABLE 1. Summary of conductance data for molecules xPhx. 
(a) Galculated for the fully-extended (all-antz conformer) 
molecule, terminated by two Au-S bonds, using Spartan 14 
molecular mechanics calculations. 


two molecules in the gap. Next, we examined 4Ph4, 
3Ph3 and IPhl. It was necessary to increase the set 
point current for these molecules, to 10, 20 and 30 nA re¬ 
spectively, since the molecules are progressively shorter, 
hence the tip has to be brought into closer proximity to 
the surface to ensure a reasonable percentage of success¬ 
ful junction formation events on tip retraction. Fig. 
shows the results obtained for 4Ph4, and Fig. those 
for IPhl; corresponding results for 3Ph3 are shown in 
Fig.3 in the SM.I^ A small proportion of the I — z 
traces from these experiments showed some evidence of a 
final, lower-conductance plateau, following the plateau 
that gave rise to the above conductance peaks. This 
suggested that lower-conductance junctions might form 
with these molecules at larger tip-substrate distances, 
but these plateaux do not give rise to a discernable sepa¬ 
rate conductance peak in either one-dimensional or two- 
dimensional histograms (see e.g. Fig. Ha) and Hb) 
where they are lost in the low-current noise at < 0.8 nS). 
These events were quite rare, even when the set-point 
current was set to a smaller value to favour the forma¬ 
tion of these longer, lower conductance plateaus over the 
shorter, higher conductance events, and so we did not 
analyse these events further. The data for all experiments 
are compared in Table |lj We used a previously-published 
procedure to measure the break-off distance (the distance 
z at which the junctions break down)P^and to correct this 
for the initial set-point distance zq (the initial height of 
the STM tip). These break-off distances were similarly 
plotted as histograms (see SM,I^Fig.6). Since, as is evi¬ 
dent from Fig. 2 and 3, the conductances of the junctions 
are rather invariant with distance for these highly flexible 
molecules, we use the 95th percentile break-off distance 
as an estimate of the maximum length of the junction 
at breakdown, and it can be seen from Table 1 that this 
correlates well with the theoretical length of a junction 
between the two contacting gold atoms, approximated for 
the fully-extended molecular conformation using molec¬ 
ular mechanics (Spartan 14). The exception is IPhl, 
where the junctions at breakdown are somewhat longer 
than expected for a molecule of this length. It is possible 
that for a very short molecule, extension of the junction 
results in gold atoms being pulled out from the surface 
prior to junction breakdown. 
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FIG. 4. Plots of In (conductance in nS) vs. molecular 
length between the contact S atoms for the xPhx series (this 
work), and for the low and medium conductance groups as 
recently measured for 1,5-pentanedithiol, 1,6-hexanedithiol, 
1,8-octanedithiol and 1,10-decanedithiol in UHV^and in am¬ 
bient (1,12-dodecanedithioP^. These examples are shown as 
these molecules have similar lengths to the shorter members 
of the xPhx series; the conductances of junctions with longer 
alkanedithiol have not been experimentally determined in this 
fashion as the conductance values would be too low. 


It is interesting to note that, alone of the molecules 
studied here, IPhl has previously been studied by the 
STM-BJ technique, and at 0.6 V bias voltage, the con¬ 
ductance measured was 80 nS.^ It has previously been 
shown that the STM-BJ method favours the formation 
of high conductance molecular junctions, and this has 
been interpreted in terms of the contact atoms binding 
to gold at step edges or similar defects, producing a junc¬ 
tion with a shorter metal — metal distance and with 
additional channels for conduction through the molecule 
due to the step edge gold atom in close proximity to the 
molecular backbone.^ In Fig.|^ we compare the conduc¬ 
tance as a function of length between the sulfur contact 
atoms (calculated for the fully-extended molecules using 
the molecular mechanics approach in Spartanl4) for the 
xPhx series with recent data for the alkanedithiols, mea¬ 
sured in UHV.I^ A notable feature of the conductance 
data for the xPhx series is the very low apparent value 
of [3 (calculated from the slopes of these lines; Eq. 0 ), 
1.52 nm 

Since the only structural change in the molecules xPhx 
is the length of the alkyl chains, it is pertinent to compare 
the (3 values with those found for alkanedithiols in related 
STM-based experiments; reported values r ange fr om 7- 
8.2 nm“^ for all three conductance groups. gi^ce 
the Fermi energy of the contacts has to fall between the 
HOMO and LUMO of the molecule in the junction, for 
molecules conducting via a tunnelling mechanism one ex¬ 
pects the p value to depend upon the HOMO-LUMO sep¬ 
aration, with more conjugated molecules having a smaller 
HOMO-LUMO gap and hence a smaller (3. For instance, 
for a series of oligophenyl methanethiols (C 6 H 5 CH 2 SH, 


C 6 H 5 C 6 H 4 CH 2 SH and C 6 H 5 C 6 H 4 C 6 H 4 CH 2 SH), p was 
determined as 4.6 ± 0.7 nm“^ using the Galn/Ga 203 
eutectic contact technique on self-assembled monolayers 
of the molecules on template-stripped gold electrodeJ^ 
(a similar value, 4.2 ± 0.7 nm“^, was determined ear¬ 
lier using a conducting AFM techniqu^^, while for 
a series of oligothiophenes with 5, 8 , 11 and 14 2,5- 
coupled thiophene rings, (3 was 1 nm“^, measured by 
the STM break junction method.l^ The lowest p val¬ 
ues so far measured have been for extended viologens 
(0.06 nm“^)P^and for conjugated fused oligoporphyrins, 
oligoporphyrin-ethynylene and oligoporphyrinbutadiyny- 
lene systems, where values of 0.2 to 0.4 nm“^ have been 
found, the exact value depending on struct ure.l^®^ It 
should be noted that although a low value of P is of¬ 
ten taken as an indication that the conductance mecha¬ 
nism could be hopping rather than superexchange (tun¬ 
nelling), it has been shown that the conductance be¬ 
haviour of even the most conjugated of these families 
of porphyrins as a function of temperature are still de¬ 
scribed well by theoretical models involving tunnelling.l^ 
In this context, a value of P oicd,. 1.5 nm“^ is remarkably 
small for the xPhx series, given that the only variable in 
this series is the number of methylene groups, for which 
we would expect p to be 7-8 nm“^. To check for a possi¬ 
ble hopping mechanism, we have carried out conductance 
measurements as a function of temperature on a repre¬ 
sentative molecule from this series, 6 Ph 6 (for data, see 
Fig.5 in the SIVP^. The conductance of this molecule 
as determined by the I{z) technique did not vary, within 
experimental limits, over the temperature range 20-100 
°C. This confirms that a hopping mechanism is not op¬ 
erative, as expected since a single benzene ring is not 
redox-active within conventional potential limits. 


III. THEORETICAL MODELLING AND 
DISCUSSION 


To gain insight into the conductance behavior of the 
xPhx system we performed DFT and electronic trans¬ 
port calculations with the help of the code GaussiaiP^ 
, using the PBE functionaP^ and the LANL2DZ basis 
set, ESI unless otherwise specified. The molecular junc¬ 
tions were built by matching the S atoms of the molecule 
with the S atoms of the S-Au structure. The elect ronic 
transport calculations were carried out with ANT.G,EiH^ 
which is built as an interface to Gaussian. This code com¬ 
putes the electronic transmission using non-equilibrium 
Green’s function techniques in the spirit of the Landauer 
formalism, employing parametrized tight-binding Bethe 
lattices in the far electrode descript ion. ESI A CRENBS 
basis set was usecP^ for the Au atoms in the outermost 
layers. Further details of the calculations involving alter¬ 
native contact sites (top, step) are given in the SM.ESl We 
optimized the molecular structures for the geometries of 
IPhl, 3Ph3, 4Ph4 and 6 Ph 6 in the gas phase and then re¬ 
laxed a S atom on top of a 19 gold atom cluster (together 
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FIG. 5. (a) Frontier orbitals of 6Ph6; the H atoms of the 
thiols are omitted as these are lost on adsorption to gold; 
(b) schematic picture of the metal-semiconductor-metal junc¬ 
tion created once the molecule is connected to the electrodes. 
The vertical arrow indicates the energy difference between 
the metal Fermi level and the valence band; (c) A model of 
the 6Ph6 molecule in a junction, forming a back-to-back pair 
of metal-semiconductor interfaces where the benzene ring is 
treated as the semiconductor bulk, the gold contacts are the 
metal and the S and outermost C atoms represent an interface 
region, (d) Schematic electronic potential along the molecular 
junction depicted in panel (c) (solid line) and along a junction 
comprising a shorter molecule (dashed line). In the case of 
the longer molecule, the potential generated by the two point 
dipoles at the interfaces gradually vanishes in the bulk region, 
which is therefore only partially affected by the presence of 
the dipoles. This is not the case for the shorter molecule, 
which does not offer enough distance for this recovery and 
where the effect of the dipoles in the bulk region is stronger. 


with the gold atoms in the top layer). Then we analyzed 
their electronic structure after removing the terminal H 
atoms, which are known to be lost upon adsorption onto 
gold under our conditions.!^ 

For all molecules, the HOMO-2 is localized on the ben¬ 
zene ring and is found to rise in energy as the chain length 
increases (see Fig. 9). As will be seen later in the dis¬ 
cussion, the length dependence of the alignment of this 
orbital with the Fermi energy plays a determining role 
in the conductance behavior. The LUMO, HOMO and 
HOMO-1 all consist of spin degenerate levels, localized 
on the S atoms. Orbitals localized on the alkyl chains are 
far below in energy. In the left panel of Fig. we show 
the orbitals for the longest molecule (6Ph6) with the cor¬ 
responding energies. It can be observed that the LUMO 
(localized on the S atoms) is quite close in energy to the 
HOMO, while the LUMO+1 (localized on the benzene) 
is far above. More importantly, a clear spatial separation 
appears between the states localized on the S atoms and 
those localized on the benzene ring. 

In our view, the LUMO+1 and HOMO-2 can be 
likened to the conduction-valence band structure in a 


semiconductor, while the states in between (LUMO, 
HOMO and HOMO-1) can be likened to the surface 
midgap states of the semiconductor as they are mainly 
localized on the S atoms that connect to the electrodes. 
Consequently, placing this molecule between two metal 
electrodes will form a junction that is expected to behave 
like a pair of back-to-back metal-semiconductor junctions 
in which the chemical terminations of the semiconductor 
create states within its band gap (as schematically de¬ 
picted in Fig. [^). For the purposes of this analogy we 
refer to the benzene core of the molecule as the “semicon¬ 
ductor bulk”. The whole junction can then be thought 
of as divided into five main blocks: metal-interface-bulk- 
interface-metal (Fig. HD- Each diode consists of a metal- 
(organic) semiconductor point-like junction. In our case, 
the semiconductor is intrinsic and the metal Fermi level 
falls between the HOMO and the LUMO of the cen¬ 
tral benzene ring (i.e. between the HOMO-2 and the 
LUMO+1 of the complete system). 

It is not straightforward, a priori, to know where the 
Fermi level will lie in this gap. However, insights can be 
gained by considering what happens at the interface be¬ 
tween a metal and a semiconductor. In analogy with the 
semiconductor surface states, the states localized on the 
S atoms of the xPhx molecules should be strongly hy¬ 
bridized when interacting with the metal, giving rise to 
the so-called ‘gateway states’.!^ The mismatch between 
the metal Fermi level and the molecular charge neutral¬ 
ity level (CNL) at the interface induces a charge transfer 
which in turn creates a point-like dipole}^ the energy lev¬ 
els in the inner part of the semiconductor bulk are then 
shifted up or down due to the dipole by a quantity which 
varies with the distance from the interface (Fig.[^). The 
role of dipole interfaces has been extensively investigated 
in the case of organic monolayers or (more generally) sin¬ 
gle molecules lying flat on surfaces^^H^ but, to the best 
of our knowledge, a thorough study concerning their roles 
in junctions comprising a single molecular nanowire held 
between two metal clusters has not been explored to the 
same extentIt is worth stressing that in this case, 
since the semiconductor is intrinsic, no depletion layer is 
formed. Consequently, variations in the energy position 
of the bulk levels must not be ascribed to the typical 
band bending associated with it in doped semiconduc¬ 
tors. Here, the only electrostatic effect is related with the 
dipole established at the interface and will be confined to 
a narrow region close to it. Notice also that while in a 
planar junction the bands are shifted at (infinitely) large 
distances,!^ in a single-point junction, such as the one 
considered here, the effect of the dipole effect is expected 
to decay rapidly along the transport direction. In what 
follows, we proceed to analyze the energetics of the four 
molecules when placed between two gold electrodes; we 
will show how studying the electrostatic balance estab¬ 
lished at the metal-molecule interface is essential for a 
thorough understanding of the level alignment, which in 
turn impacts on charge transport across the junction. 

We constructed junctions in which each molecule is 
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FIG. 6. Representative dithiols studied theoretically, show¬ 
ing different models for adsorption sites on gold, (a): do- 
decanedithiol bound at hollow sites, (b) IPhl bound at step 
edge sites, (c) IPhl bound at hollow sites and (d) IPhl bound 
at top sites. 


connected to two gold clusters, in the hollow, top or step 
binding geometry (see Fig. |^; however, in what follows, 
we will restrict discussion to the hollow position. For 
comparison, we also studied a series of alkanedithiols of 
four different lengths (C2, C6, C8 and C12), also bound 
in the hollow position. Although their total lengths are 
different from the xPhx molecules, this comparison was 
made to give insight into the presence of any additional 
state(s) and into differences in the length dependence. 
In Fig. we show the low-bias transmission curves ob¬ 
tained for the Cx and xPhx series in the hollow geome¬ 
try. For the xPhx molecules, the presence of additional 
peaks is evident. These Breit-Wigner resonances ap¬ 
proach the Fermi level as the molecul ar len gth increases 
(consistent with previous calculation^^^^^ and they do 
not appear in the case of the Cx molecules. The pres¬ 
ence of these states creates an indentation in the barrier, 
confirming our earlier suggestion.!!^ A similar effect has 
previously been shown in another system.!^ These res¬ 
onances stem from the HOMO-2, which is localized on 
the benzene ring, thus they become narrower as the in¬ 
creasing benzene-gold distance decouples them from the 
gold contacts. To some extent, the benzene ring behaves 
as a quantum dot, weakly coupled to the electrodes by 
the alkyl chains. The trend in energy of the peaks ap¬ 
proximately follows that in the gas phase for 3Ph3, 4Ph4, 
and 6Ph6, whereas the corresponding peak of IPhl lines 
up further down in energy (at -1.48 eV). Indeed, in the 
case of this shorter molecule, exceptional behavior is ex¬ 
pected due to the strong hybridization of the orbitals. In 
the case of the Cx molecules, the states closest to the 
Fermi level are localized mainly on the S atoms, while 
states localized on the whole chain appear only below -2 
eV. 

To discern the different contributions in the transmis¬ 
sion curves of the xPhx molecules, we calculated the DOS 
projected onto the molecular orbitals. The results are il¬ 



FIG. 7. Transmission as a function of energy for Gx and 
xPhx in the hollow binding geometry and DOS projected onto 
HOMO -2 and HOMO-3 for 6Ph6 in the hollow geometry. 


lustrated for 6Ph6 in the lowest panel of Fig. These 
projections clearly show that the sharp peak visible in 
the transmission curve close to the Fermi level for 6Ph6 
is due to the HOMO-2, while the HOMO-3, also mainly 
localised on the benzene ring, does not yield any peak in 
the transmission curve apart from a very small Fano res¬ 
onance (interference effects in transport through benzene 
rings have been widely studied^. The broad bumps at 
around -1 eV derive from the gateway states, i.e. from 
the hybridization of S-Au states at the interface. Notice 
that using the same voltage 0.6 V as in the experiments 
is not expected to change the overall picture. To con¬ 
firm this, we show, in the SM, a comparison between the 
transmission curves for 4Ph4 calculated at 0 and 0.6 V: 
the position of the HOMO-2 is not affected, while the bias 
voltage seems to rather have an influence on the interface 
states. 

We now turn to analyse these systems quantitatively. 
In what follows, we will focus on the metal-molecule in¬ 
terface of our system in order to discern the various con¬ 
tributions that determine the level alignment in the stud¬ 
ied junctions. We again focus on the hollow binding ge¬ 
ometry. To begin with, we have first analyzed the Mul- 
liken charges in order to discern whether any significant 
charge transfer takes place in this system (see Table I 
in the SIVP^. They show an increased electron charge 
(around 0.3) on each gold cluster and corresponding pos¬ 
itive charge on the outermost carbon atoms, indicating 
charge transfer from the molecule onto the metal. Sub¬ 
sequently, we have analysed a system consisting of each 
molecule connected to a single Aul9 cluster in a hollow 
position, with the aim of understanding the level align¬ 
ment at just one metal-molecule interface. In this case. 
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as well as the standard DFT calculation, we have also 
performed constrained DFT (c-DFT) calculations, as ex¬ 
plained below. By connecting a Bethe lattice to the metal 
cluster, a new system is built, where the molecule is now 
connected to a semi-infinite electrode (see Ref. [42] for 
more details). The system can be hence divided into two 
parts, namely the lead and the main region that com¬ 
prises the molecule and a certain number of metal layers. 
The density matrix of this region is calculated as 

2 7^ 

P(/i) =- Im dPG+(P;/i) (2) 

^ J—oo 

where the retarded Green’s function is given by 

G+{E- ii) = [{E-ii)-S-H- E(E, (3) 

Here, H is the Hamiltonian of the main region, S is the 
lead self energy, which describes the coupling to the semi¬ 
infinite lead, and /i is a quantity by which the on-site 
energies of the Hamiltonian must be shifted in order to 
ensure the total-charge neutrality (/i is opposite in sign to 
the chemical potential). In our c-DFT any charge trans¬ 
fer between the two parts is forbidden. Both surface and 
molecule were forced to keep a number of electrons as 
to maintain charge neutrality. To this aim, /i was com¬ 
puted in alternation so as to meet the imposed charge 
constraints on the metal and the molecule. If Ns and 
Nmol are the desired number of electrons in the surface 
and the molecule, respectively, then two different chemi¬ 
cal potentials fis and fimoi are calculated so that 

NAOs 

Tr[Ps{^ls)Ss\= Y. P^^G)S{i,3)=Ns (4) 

7i=i 

N AO s+N AOjnoi 

7^'^[-PmoZ (/^moz)^moz] ~ ^ ^ j)^(b4) ~ ^moh 

i,j=NAOs + l 

( 5 ) 

where NAOs and N AO mol are the number of atomic 
orbitals of the surface and the molecule, respectively. At 
each iteration step the density matrix is then built as 
a block matrix as follows: 

po = ( d ( 6 ) 

VPS'-moZ ^moiy^mol) J 

The off-diagonal terms in the submatrices Pg-moi were 
set to zero. The so-built density matrix is out of equilib¬ 
rium, with two different chemical potentials in the metal 
and the molecule (at approximately -5.96 eV and -4.2 eV 
for metal and molecule, respectively), consequently giv¬ 
ing rise to a step potential across the system. The Fermi 
level of the molecule calculated in this way can be inter¬ 
preted as the charg e neutrality level (CNL) discussed in 
the hterature,^^^^ but obtained in a more rigorous way. 

It defines the direction in which electronic charge is most 



FIG. 8. On-site energies (averaged over the atomic orbitals of 
each atom) in the three innermost An atoms and the molecu¬ 
lar C atoms in the Aul9-xPhx systems for all four molecules, 
in both constrained and unconstrained DFT. 


likely to be transferred. In our case, the GNL is definitely 
higher than the An Fermi level, so charge transfer from 
the molecule onto the metal is expected, confirming what 
was indicated by the Mulliken charges in the junction. 

We now turn to analyze what happens once the elec¬ 
tronic charge has been transferred. To this end, we have 
computed the on-site energies along the Aul9-molecule 
junctions and compared them with those resulting from 
the constrained calculation (the results are shown in 
Fig. 1^. They were evaluated as the average (over the 
atomic orbitals of each atom) of the diagonal terms in 
the Kohn-Sham Hamiltonian. Since the absolute values 
along the vertical axes are not physically relevant, as they 
are strongly related to the basis-set choice for each atom, 
they have been offset for clarity. Hence, we will only fo¬ 
cus on the comparison between the values obtained in 
two kinds of calculations. At the interface, we see a clear 
sign of the expected dipole; in the unconstrained calcu¬ 
lation the values are shifted upwards in the metal and 
downwards in the molecule. In the latter, the effect is 
stronger in the carbon atoms closer to the metal, fading 
away as the distance from the metal increases. In IPhl, 
the shortest molecule, the benzene ring is quite close to 
the gold and shows some positive charge (0.1) in the un¬ 
constrained calculation, as compared to that constrained. 
Indeed, for IPhl the on-site energies are strongly affected 
also in the G atoms of the benzene ring. It is also worth 
adding that, when a second cluster is added at the other 
end of the molecules to complete the junctions, the effects 
caused by the dipole in the bulk region are expected to 
increase. The Mulliken charges along the Aul9-molecule 
junctions in the constrained and unconstrained calcula¬ 
tions are shown in the SM.I^ 

We summarize all these results as follows: when the 
xPhx molecules approach the metal, Fermi-level align¬ 
ment is achieved upon charge transfer from the sulfur, or 
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Molecular length between contact S atoms (nm) 


FIG. 9. Energetic position of the HOMO-2 for each molecule 
in the junction and in the gas phase. All values have been 
shifted to set those for 6Ph6 to zero. It can be seen that 
the original displacement of the HOMO-2 in the gas phase 
increases in the junction. 


FIG. 10. The calculated In(conductance) as a function of 
the molecular length for Gx and xPhx in the hollow binding 
geometry 


the outermost carbon atoms, onto the metal. This gen¬ 
erates a dipole potential mainly near the interface. The 
molecular bulk states (HOMO-2 levels) then behave like 
spectators, with their energetic positions being dragged 
down under the influence of the dipole potential accord¬ 
ing to the distance from the interface (Fig. ^), and never 
crossing the Fermi level (see, for example, the case of 
the top binding geometry for the hypothetical 8Ph8 in 
the SlVpl). The level alignment of the bulk states is 
governed purely by electrostatics because of the clear 
separation from the interface (apart from IPhl) as al¬ 
ready discussed. Note that this would not be possible 
in fully conjugated molecules, where other effects such 
as hybridization would also play a role. Notice that 
any electrostatic-induced shift is added to the one the 
HOMO-2 already presents in the gas phase for increas¬ 
ing chain length (see, in Fig.[^ a comparison between the 
energy shift experienced by this orbital in the gas phase 
and in the junction). 

It should be recalled that there is some uncertainty 
in the position of the HOMO-2 peak in the transmis¬ 
sion curves as it corresponds to a Kohn-Sham state, and 
it might be located at lower energies. Nevertheless, the 
computed transmission curves provide an insight into the 
alignment dynamics, i.e. they clearly show how the align¬ 
ment of the bulk states is governed by the gas-phase be¬ 
havior, modulated by what happens at the interface. To 
support our interpretation of the numerical results, it is 
important to verify that the benzene ring falls inside the 
space region in which the dipole-potential effect is con¬ 
fined. To this aim, we have calculated the potential gen¬ 
erated in the middle of the benzene ring by the charges 
at the interface. The results are reported in Table II 
of SlvP and show decreasing values for increasing chain 
length, with maximum values of around 1 eV. 

Overall, the scenario described above provides a good 


rationalisation for the experimentally observed low value 
of /3. There are two counteracting effects, namely the ex¬ 
pected decay due to increasing molecular length on the 
one hand, together with the increasing sharpness of the 
HOMO-2 resonance in the molecules with longer poly¬ 
methylene bridges, and the approach to the Fermi energy 
of the resonance peak on the other hand. Fig. pM| shows 
a plot of calculated zero bias conductance vs. molec¬ 
ular length for the xPhx and Cx molecules. Although 
the absolute values of the conductances differ substan¬ 
tially from the experimentally observed values, it can be 
clearly seen that the calculated (3 value is smaller for the 
xPhx series than for the Cx molecules, although the fac¬ 
tor by which they differ (ca. 1.5) is substantially smaller 
than that seen experimentally (ca. 4.8). Nevertheless, 
given the known problems with DFT, and the limita¬ 
tions inherent in the modeling of the contacts, we believe 
that this offers a good explanation for the experimental 
observations. Ultimately, possible asymmetries in the ex¬ 
perimental junctions should be taken into account (not 
considered here) which lead to transmission resonances 
with a maximum height lower than 1. While this is not 
expected to change the conductance values for the Cx 
molecules, for which the resonances show up further be¬ 
low the Fermi energy, it might change the conductance 
values for the xPhx molecules and, consequently, the {3 
factor. 


IV. CONCLUSIONS 

We studied the series of molecules 
HS(CH2)^[1,4C6H4](CH2)^SH (x = 1, 3, 4, 6) and 
found that the rate of conductance decay with molecular 
length is smaller than that measured previously for alka- 
nedithiols. This was found to be due to the decreasing 
distance in energy of the Breit-Wigner resonance from 
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the Fermi level for increasing molecular length. This 
trend originates from the orbital behaviour in the gas 
phase but is modulated by the electrostatic balance 
established at the electrode contact region. 


V. ACKNOWLEDGMENTS 

This research was supported by the EPSRC (Grant 
No. EP/H035184/1), by MINECO under Grant No. 
EIS2013-47328, by the European Union structural funds 
and the Comunidad de Madrid MAD2D-CM Program 
under Grant. P2013/MIT-2850, and by Generalitat Va- 
lenciana under Grant PROMETEO/2012/011. We thank 
M. Soriano for help with the computational part and ac¬ 
knowledge the computer resources, technical expertise, 
and assistance provided by the Centro de Computacion 
Cientifica of the Universidad Autonoma de Madrid. 


^ R. J. Nichols, W. Haiss, S. J. Higgins, E. Leary, S. Martin, 
and D. Bethell, Physical Chemistry Chemical Physics 12, 
2801 (2010). 

^ F. Chen, J. Hihath, Z. Huang, X. Li, and N. Tao, Anna. 

Rev. Phys. Chem. 58 , 535 (2007). 

^ M. S. Hybertsen, L. Venkataraman, J. E. Klare, A. C. 
Whalley, M. L. Steigerwald, and C. Nuckolls, Journal of 
Physics: Condensed Matter 20 , 374115 (2008). 

^ F. Chen and N. Tao, Accounts of Chemical Research 42 , 
429 (2009). 

^ J. C. Cuevas and E. Scheer, Molecular electronics: an 
introduction to theory and experiment (World Scien¬ 
tific,Singapore, 2010). 

® J. Park, A. N. Pasupathy, J. 1. Goldsmith, C. Chang, 
Y. Yaish, J. R. Petta, M. Rinkoski, J. P. Sethna, H. D. 
Abruha, P. L. McEuen, et a/.. Nature 417, 722 (2002). 

^ J. Reichert, R. Ochs, D. Beckmann, H. Weber, M. Mayor, 
and H. von Lohneysen, Phys. Rev. Lett 88, 176804 (2002). 
® B. Xu and N. J. Tao, Science 301 , 1221 (2003). 

^ S. Y. Quek, M. Kamenetska, M. L. Steigerwald, H. J. 
Choi, S. G. Louie, M. S. Hybertsen, J. Neaton, and 

L. Venkataraman, Nature Nanotechnology 4, 230 (2009). 

M. Kamenetska, S. Y. Quek, A. Whalley, M. Steiger¬ 
wald, H. Choi, S. G. Louie, C. Nuckolls, M. Hybertsen, 
J. Neaton, and L. Venkataraman, Journal of the Ameri¬ 
can Chemical Society 132 , 6817 (2010). 

F. Chen, X. Li, J. Hihath, Z. Huang, and N. Tao, Journal 
of the American Chemical Society 128 , 15874 (2006). 

L. Venkataraman, J. E. Klare, 1. W. Tam, C. Nuckolls, 

M. S. Hybertsen, and M. L. Steigerwald, Nano Letters 6, 
458 (2006). 

S. Martin, W. Haiss, S. Higgins, P. Cea, M. C. Lopez, and 
R. J. Nichols, The Journal of Physical Chemistry C 112, 
3941 (2008). 

Y. S. Park, A. C. Whalley, M. Kamenetska, M. L. Steiger¬ 
wald, M. S. Hybertsen, C. Nuckolls, and L. Venkataraman, 
Journal of the American Chemical Society 129 , 15768 
(2007). 

R. Parameswaran, J. Widawsky, H. Vazquez, Y. Park, 
B. Boardman, C. Nuckolls, M. Steigerwald, M. Hybertsen, 
and L. Venkataraman, The Journal of Physical Chemistry 
Letters 1 , 2114 (2010). 

W. Haiss, H. van Zalinge, S. J. Higgins, D. Bethell, 
H. Hobenreich, D. J. Schiffrin, and R. J. Nichols, Journal 
of the American Chemical Society 125 , 15294 (2003). 


W. Haiss, S. Martin, E. Leary, H. v. Zalinge, S. J. Higgins, 
L. Boufiier, and R. J. Nichols, The Journal of Physical 
Chemistry C 113 , 5823 (2009). 

E. Leary, S. J. Higgins, H. van Zalinge, W. Haiss, and 
R. J. Nichols, Chem. Commun. , 3939 (2007). 

G. K. Ramachandran, J. K. Tomfohr, J. Li, O. F. Sankey, 

X. Zarate, A. Primak, Y. Terazono, T. A. Moore, A. L. 
Moore, D. Gust, et a/.. The Journal of Physical Chemistry 
B 107 , 6162 (2003). 

B. Q. Xu, X. L. Li, X. Y. Xiao, H. Sakaguchi, and N. J. 
Tao, Nano Letters 5 , 1491 (2005). 

R. Yamada, H. Kumazawa, T. Noutoshi, S. Tanaka, and 

H. Tada, Nano letters 8, 1237 (2008). 

S. K. Lee, R. Yamada, S. Tanaka, G. S. Chang, Y. Asai, 
and H. Tada, ACS Nano 6, 5078 (2012). 

E. J. Dell, B. Capozzi, K. H. DuBay, T. C. Berkelbach, 
J. R. Moreno, D. R. Reichman, L. Venkataraman, and 

L. M. Campos, Journal of the American Chemical Society 
135 , 11724 (2013). 

G. Sedghi, L. J. Esdaile, H. L. Anderson, S. Martin, 

D. Bethell, S. J. Higgins, and R. J. Nichols, Advanced 
Materials 24 , 653 (2012). 

G. Sedghi, V. M. Garcfa-Suarez, L. J. Esdaile, H. L. An¬ 
derson, C. J. Lambert, S. Martin, D. Bethell, S. J. Higgins, 

M. Elliott, N. Bennett, et al., Nature Nanotechnology 6, 
517 (2011). 

G. Sedghi, K. Sawada, L. J. Esdaile, M. Hoffmann, H. L. 
Anderson, D. Bethell, W. Haiss, S. J. Higgins, and R. J. 
Nichols, Journal of the American Chemical Society 130 , 
8582 (2008). 

Z. Li, T.-H. Park, J. Rawson, M. J. Therien, and 

E. Borguet, Nano Letters 12, 2722 (2012). 

C. Wang, A. S. Batsanov, M. R. Bryce, S. Martn, R. J. 
Nichols, S. J. Higgins, V. M. Garca-Surez, and C. J. Lam¬ 
bert, Journal of the American Chemical Society 131 , 15647 
(2009). 

“See supplemental material at [url will be inserted by 
publisher] for details about synthetic procedures, single 
molecule experiments, theoretical calculations and spectra 
of compounds,”. 

W. Haiss, D. Lackey, J. K. Sass, and K. H. Besocke, The 
Journal of Chemical Physics 95 , 2193 (1991) 

X. Xiao, B. Xu, and N. J. Tao, Nano Letters 4, 267 (2004). 
E. Pires, J. E. Macdonald, and M. Elliott, Nanoscale 5, 
9397 (2013). 



10 


Z. Li, L Pobelov, B. Han, T. Wandlowski, A. Blaszczyk, 
and M. Mayor, Nanotechnology 18 , 044018 (2007). 

A. Busiakiewicz, S. Karthauser, M. Homberger, 
P. Kowalzik, R. Waser, and U. Simon, Physical 
Chemistry Chemical Physics 12, 10518 (2010). 

D. Fracasso, M. I. Muglali, M. Rohwerder, A. Terfort, and 
R. C. Chiechi, The Journal of Physical Chemistry C 117 , 
11367 (2013). 

D. J. Wold, R. Haag, M. A. Rampi, and C. D. Frisbie, 
The Journal of Physical Chemistry B 106 , 2813 (2002). 

V. Kolivoska, M. Valasek, M. Gal, R. Sokolova, J. Bulick- 
ova, L. Pospisil, G. Meszaros, and M. Hromadova, The 
Journal of Physical Chemistry Letters 4, 589 (2013). 

M. J. Frish, G. W. Trucks, and H. B. S. et ah, GAUSSIAN 
03, Revision B.01,Gaussian. Inc., Pittsburg (2003). 

J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. 
Lett 78 , 1396 (1997). 

W. R. Wadt and P. J. Hay, The Journal of Chemical 
Physics 82 , 284 (1985). 

J. J. Palacios, D. Jacob, A. J. Perez-Jimenez, E. S. Eabian, 

E. Louis, and J. A. Verges, ALACANT ab initio quantum 
transport package, see URL: http://alacant.dfa.ua.es. 

D. Jacob and J. J. Palacios, Journ. Chem. Phys. 134, 
044118 (2011). 

R. Ross, J. Powers, T. Atashroo, W. Ermler, L. LaJohn, 
and P. Christiansen, The Journal of chemical physics 93 , 
6654 (1990). 

J. C. Love, L. A. Estroff, J. K. Kriebel, R. G. Nuzzo, and 
G. M. Whitesides, Chemical Reviews 105 , 1103 (2005). 

J. R. Widawsky, W. Chen, H. Vzquez, T. Kim, R. Breslow, 
M. S. Hybertsen, and L. Venkataraman, Nano Letters 13 , 


2889 (2013). 

E. Elores, J. Ortega, and H. Vazquez, in Advances in Solid 
State Physics (Springer, 2009) pp. 207-218. 

G. Heimel, L. Romaner, J.-L. Brdas, and E. Zojer, Surface 
Science 600 , 4548 (2006) 

M. Bokdam, D. Qakir, and G. Brocks, Applied Physics 
Letters 98, 113303 (2011). 

O. T. Hofmann, D. A. Egger, and E. Zojer, Nano letters 
10 , 4369 (2010). 

S. Braun, W. R. Salaneck, and M. Eahlman, Advanced 
Materials 21, 1450 (2009). 

R. Gutierrez, G. Eagas, K. Richter, E. Grossmann, and 
R. Schmidt, EPL (Europhysics Letters) 62 , 90 (2003). 

Y. Xue, S. Datta, and M. A. Ratner, The Journal of Chem¬ 
ical Physics 115 , 4292 (2001). 

G. Peng, M. Strange, K. S. Thygesen, and M. Mavrikakis, 
The Journal of Physical Chemistry C 113, 20967 (2009). 
E. Leonard and J. Tersoff, Physical Review Letters 84, 
4693 (2000). 

E. Hiiser, Diploma Thesis, Universitat Karlsruhe (2009). 

C. Benesch, M. Cizek, J. Klimes, 1. Kondov, M. Thoss, 
and W. Domcke, The Journal of Physical Chemistry C 
112, 9880 (2008). 

L. A. Zotti, E. Leary, M. Soriano, J. C. Cuevas, and J. J. 
Palacios, Journal of the American Chemical Society 135 , 
2052 (2013). 

D. Darau, G. Begemann, A. Donarini, and M. Grifoni, 
Physical Review B 79 , 235404 (2009). 

H. Vzquez, E. Elores, and A. Kahn, IPAP Conf. Series 6 
, 1 (2005). 


